arXiv: 1505.00216v4 [astro-ph.GA] 21 Feb 2016 


Draft version February 23, 2016 

Preprint typeset using DTgX style emulateapj v. 5/2/11 


X-RAY SOURCES IN THE DWARF SPHEROIDAL GALAXY DRACO 

E. Sonbas 1,2 , B. Rangelov 2 , O. Kargaltsev 2 , K.S. Dhuga 2 , J. Hare 2 , and I. Volkov 2,3 
1 University of Adiyaman, Department of Physics, 02040 Adiyaman, Turkey 
department of Physics, The George Washington University, Washington, DC 20052, USA and 
3 College of Computer, Mathematical, and Natural Sciences, University of Maryland, College Park, MD 20742, USA 

Draft version February 23, 2016 

ABSTRACT 

We present the spectral analysis of an 87 ks XMM-Newton observation of Draco, a nearby dwarf 
spheroidal galaxy. Of the approximately 35 robust X-ray source detections, we focus our attention on 
the brightest of these sources, for which we report X-ray and multiwavelength parameters. While most 
of the sources exhibit properties consistent with AGN, few of them possess characteristics of LMXBs 
and CVs. Our analysis puts constraints on population of X-ray sources with Lx > 3 x 10 33 erg s -1 in 
Draco suggesting that there are no actively accreting BH and NS binaries. However, we find 4 sources 
that could be LMXBs/CVs in quiescent state associated with Draco. We also place constraints on 
the central black hole luminosity and on a dark matter decay signal around 3.5 keV. 

Subject headings: galaxies: individual (Draco dwarf) — galaxies: dwarf — X-ray: binaries — novae, 
cataclysmic variables — stars: carbon — dark matter 


1. INTRODUCTION 

The large spirals of the Local Group, the Milky Way 
and Andromeda (M31), are known to have a significant 
number of faint satellite galaxies; some of the faintest 


vey (SPSS; e.g., 


Newberg et al. 

2002 

Willman et al. 

spheroidals (dSp 

ns) are 

by far the most 


numerous class among these satellites and typically the 
least luminous. In general, these objects are gas-poor, 
spatially diffuse, and exhibit signific antly elevated mass- 
to-light ratios (M /L ~ few x 100; |Mateo||1998| |Baum- 
gardt et al. |2008|) . The latter implies that the dSphs 
occupy an important position in the mass spectrum rep¬ 
resenting the smallest known structures where dark mat¬ 
ter dominates. The dSphs have low mean metallicity, and 
thus their study can provide a probe of stellar evolution 
at these low metallicities. In addition, these objects are 
of interest because they may have abundance rati os t hat 
are different from those of the Milky Way (Kirby et al. 
2011a[bj) 


A recent study on the Sculptor dwarf galaxy reported 
the detec tion of five X-ray binari es with Lx > 6 x 10 33 


erg s " (jMaccarone et al. 2005). The number seems 


to be surprisingly high given the expectations based 
purely on scaling and integrated properties such as 
central density, core radius, and velocity dispersion. 
However, this number appears to be in agreement with 
predictions of the model proposed by Piro et al. (2002), 
where the discovered X-ray sources are potentially low 
duty cycle transients. By assuming typical velocities for 
LMXBs, Dehnen & King (2006) used the X-ray data to 
place a mass limit of > 10 9 M 0 on the amount of dark 
matter (DM) in Sculptor required to retain the LMXBs. 
In addition, the authors note that there should be an 
extended halo of qLMXBs which may be observable. 
Simulation studies, involving tidal stripping, (Read et 
al. 2006b), suggest comparable properties for Draco. 
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Located at ~ 80 kpc, Draco is a faint, metal-poor 



equilibrium, the high stellar velocity dispersion implies 
an extremely high mass-to-light ratio (M/L) of about 
146 ±42. In order to probe whether the main drivers 
of scaling in dSphs are the internal collisions and en¬ 
counters, which strongly depend on the (core) density of 
these systems, a population census of LMXBs, qLMXBs 
and cataclysmic variables (CVs) is highly desirable. 
This was one of the main goals of our XMM-Newton 
observation whose results we report in this paper. 

The paper is organized as follows. Section 2 reports 
the details of observations, the data reduction and the 
source selections criteria. In Section 3, we outline our 
source classification methodology. In Section 4, the spec¬ 
tral analysis is described for the interesting sources, lim¬ 
its on the X-ray emission from the central black hole 
are presented, and the diffuse emission spectrum near 
3.5 keV is presented. In section 5, we summarize our 
results. 


2. OBSERVATIONS AND DATA ANALYSIS 
2.1. XMM-Newton 

The 87 ks image of Draco (Fig. 1; RA = 17 /l 20 m 12.4 s 
and DEC = ±57°54'55 / /3) was obtained by the Euro¬ 
pean Space Agency’s (ESA) X-ray Multi-Mirror Mission 
- Newton ( XMM-Newton ) between August 4 - 28, 2009 
(PI: K. Dhuga; see Table 1 for details). 

We used the 3XMM-DRf0 (3XMM) catalog (Rosen 


et al. 2015) for the automated multi wavelength (MW) 
classification described in Section 3.1. The parameters 
obtained from the 3XMM catalog are the source coordi¬ 
nates and X-ray fluxes in 4 different energy bands. 


1 http://xmmssc.irap.omp.eu/Catalogue/3XMM-DR5/3XMM- 
DR5_Catalogue_User_Guide.html±DocHistory 
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Fig. 1. — Left: EPIC image of Draco smoothed with gaussia m r = 2" kernel. The HST pointings are shown in wh ite. With a tidal rad ius 
of ~ 40 ; , Draco extends beyond the XMM field of view (lOdenkirchen et al.||200l|). The half-light radius of ~10 / (jGilmore et al.||2007| is 
shown with white circle. Right: Color-coded (Blue: u , Green: g , and tied: r) bDSS image of Draco. The X-ray sources are shown with 
green crosses, and numbered according to Table 3. Both images show the same region of the sky (north is up, east is to the left). 



Fig. 2. — HST cutouts (8 ;/ X 8") centered on the X-ray sources observed by HST. Each circle has a 2" radius, which represents the 
typical positional uncertainty for these XMM-Newton sources. Note that for source 1, a very bright foreground star is also seen in the HST 
images, and is not shown because the image would be flooded with the light from the star. 


For several of the brighter X-ray sources and for the 
diffuse background emission (see Section 4) we performed 
spectral analyses and extracted spectra from the EPIC 
data. The observation data files (ODFs) were processed 
using standard procedures of the XMM-Newton Science 
Analysis System (XMM-SAS version 13.0.1). The stan¬ 
dard SAS tool edetect-chain was used to perform 
source detections. We extracted spectra using circular 
(r = 17") regions centered on the X-ray sources, for each 
observation. Background subtraction was performed us¬ 
ing source-free apertures located nearby. A standard 
event screening described in the multixmmselect man- 
uaQwas applied. The spectral analysis of the X-ray data 
was performed using XSPEC version 12.7. The pn and 
MOS spectral energy channels were grouped to have at 
least 10 counts per bin. Each spectrum was fitted in the 
0.2—10 keV energy range. 


2 http://xmm.esac.esa.int/sas/current/documentation/threads/ 
multiexposures.thread.shtmi 


TABLE 1 

X-RAY OBSERVATIONS OF DRACO 


Date 

Obsld 

Observatory 

Exp a 

2009 Aug 04 

0603190101 

XMM-Newton 

17 

2009 Aug 06 

0603190201 

XMM-Newton 

18 

2009 Aug 08 

0603190301 

XMM-Newton 

16 

2009 Aug 20 

0603190401 

XMM-Newton 

18 

2009 Aug 28 

0603190501 

XMM-Newton 

18 


a Exposure in ks. 


2.2. Hubble Space Telescope 

The Hubble Space Telescope ( HST ), partly observed 
Draco with the WFPC2, ACS/WFC or WFC3/UVIS 
cameras (Table 2 lists the images used for this study). 
Figure 1 shows the HST pointings, in white, overlaid 
on top of the EPIC image. Figure 2 shows 8" x 8" 
cutouts centered on X-ray sources observed by HST. We 
use the high-resolution HST images to distinguish be¬ 
tween AGNs and stars. 
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Fig. 3. — HR diagram where sources from the Draco field are 
shown by asterisks with numbers (correspond to those in Table 3). 
The population of AGNs from our training dataset are shown in 
grey, Stars in yellow, YSOs in orange, LMXBs in cyan, HMXBs in 
green, CVs in red, and pulsars in blue. The average uncertainty of 
the sources in Table 3 is shown in the top left corner. 


TABLE 2 

HST OBSERVATIONS OF DRACO 


Date 

Instrument 

Filter 

Obsld 

Exp. a 

Res. b 

2001 Aug 18 

WFPC2 

F555W 

9043 

680 

0.1 

2004 Oct 31 

ACS C 

F606W 

10229 

8170 

0.5 

2004 Oct 30 

ACS C 

F555W 

10229 

8170 

0.5 

2004 Oct 29 

ACS C 

F606W 

10812 

7200 

0.5 

2012 Oct 26 

WFC3 d 

F606W 

12966 

3152 

0.46 

2013 Oct 14 

WFC3 d 

F606W 

12966 

3076 

0.46 


a Exposure time in seconds. 
b Pixel scale in "/pixel. 
c WFC camera. 
d UVIS camera. 


3. MULTIWAVELENGTH ANALYSIS 
3.1. Classification 

Traditional multi wavelength (MW) classification in¬ 
cludes heuristic examination of the MW source proper¬ 
ties, such as X-ray hardness ratios, X-ray/optical/NIR 
flux ratios, color-color, and color-magnitude diagrams in 
the optical and IR such as those shown in Figures 3, 4 and 
5 respectively. These diagrams allow one to devise crude 
criteria to discriminate between the sources of different 


nature (e.g., Kaplan et al. 2006| 

Misanovic et al. 

2010; 

Kargaltsev et al. 2012; Lin et al. 

2012). WhWtFeMW 

analysis is often the only reliable way to unveil the na- 


ture of unknown X-ray sources, traditional classification, 
based on drawing dividing lines in the MW diagrams is 
a very laborious process lacking any quantitative classi¬ 
fication confidence criterium. 

There are over a hundred X-ray sources in the 3XMM 
catalog in the Draco field. We visually inspected all 
sources, and after excluding source duplications (in 5 
separate observations of Draco) and spurious detections 
(e.g., edge of the chip, bad CCD columns) we selected 
35 bright X-ray sources with S/N > 6 for further inves- 

3 http://xmmssc-www.star.le.ac.uk/Catalogue/3XMM-DR4/ 
UserGuide_xmmcat.html 


tigation. Table 3 lists the source parameters: X-ray flux 
^ 0 . 5—12 in 0.5—12 keV range and two hardness ratios HR2 
and HR4 defined as HR2 = (F( 1i2 _ 2 -F 0 . 2 -i. 2 )/(^i. 2-2 + 

^o. 2 -i. 2 ) and HRA = (F 2 _ 7 — F 0 . 5 — 2 )/(^ 2 —r + ^ 0 . 5 - 2 ), 
where F x _ y are the observed fluxes in the respective en¬ 
ergy bands. Figure 3 shows the hardness ratio diagram 
for the sources in Draco compared to different classes of 
literature verified sources used for automated classifica¬ 
tion (see A1 for details). The two chosen hardness ratios 
allow for a some separation between AGNs and the rest 
of the sources, although there is an expected overlap be¬ 
tween the AGNs, CVs, and pulsars. 

Similarly, we used color-color diagrams to compare the 
sources in Draco to known sources of different classes 
(AGN, CV, pulsars etc.). The left panel of Figure 4 shows 
a good degree of separation between the stars and AGN 
in terms of IR to X-ray flux raticj^J which is often used 
to separate the two classes. The color-color diagram in 
the right panel of Figure 4 also shows clear separation 
between the AGN and stars but it also allows to see no¬ 
ticeable separation between the AGN and CVs. We also 
utilize color-magnitude diagrams (CMDs; see Figure 5) 
and visually inspect the HST and SDSS images. The 
former allows us to estimate the age of stellar sources, 
which in turn helps to determine whether a given star 
is part of Draco. If a star, whose absolute magnitude is 
calculated for the Draco distance, happens to land onto 
a < 1 Gyr isochrone, it is unlikely to be part of Draco, 
because the stellar popu latio n of th e dSph is at least sev¬ 
eral billion years old (Segall et al. 2007). AGNs have a 
spectral energy distribution (BED) that is different from 
that of stars, and, therefore, on different CMDs, AGNs 
“jump” between different isochrones, or do not land on 
one at all. 

Since we are interested in objects that belong to Draco 
we attempt to eliminate the Galactic stars. One helpful 
constraint is the optical brightness of these objects at 
the known distance to Draco. The distance modulus of 
Draco (80kpc) is ^ 19.5. Figure 5 clearly shows that we 
do not expect to see any stars brighter than M g « — 1 (a 
very conservative estimate, which translates to apparent 
magnitude m g « 18), assuming a few Gyr as the age of 


Draco. This means that stars brighter than m g 


18 


likely belong to our Galaxy rather than Draco. (This 
criteria does not apply to QSOs/AGNs, which can be 
brighter.) Sources #1, #2, #4, #11, #31 and #33 ap¬ 
pear to coincide with stars associated with our Galaxy. 
Inspection of the HST images confirms that source #1 
also contains a very bright foreground star. 

In addition, we utilize the superb angular resolution 
of HST to identify AGNs (we also check the lower- 
resolution SDSS images when HST image is not avail¬ 
able). Unfortunately, only a very limited (often single¬ 
band) HST coverage of Draco currently exists (see Fig¬ 
ure 1, left panel). Nonetheless, high resolution HST im¬ 
ages allow us to eliminate several foreground and back¬ 
ground sources to Draco (see Table 3). 

Sources #5 and #6 appear to be point-like in the HST 
images, but given their optical magnitudes (of 19-20) we 
cannot exclude the possibility of QSOs (which could also 


4 Note that in Figure 4 (left panel) we plot the spectral flux 
ratios which differ from the flux ratios by the factor equal to the 
inverse ratio of the corresponding mid-band frequencies. 
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log(F 0 / F W1 ) W2-W1 

Fig. 4. — The flux ratio and color-color diagrams traditionally used for source classification. These also represent of the many possible 
slices in MW parameter space used by the automated algorithm to create the decision tree (see Appendix). Left: X-ray-to-IR spectral 
flux (at mid-band frequency) ratio, To.5_8kev/-Gvi> vs. optical-to-IR spectral flux ratio, F g /Fw 1 , is shown in logarithmic scale. Right: 
Optical color, g — i vs. IR color, W2 — Wl. On both panels the population of AGNs from the training dataset (see Appendix) is shown in 
grey, YSOs in orange, stars in yellow, and CVs in red. The sources from the Draco field (only those with the optical and NIR counterparts) 
are shown as black asterisks and numbered according to Table 3. 



G-R 


Fig. 5. — Color-magnitude diagram showing absolute magnitude 
Mq vs. color G-R. Isochrones are plotted at 300 Myr, 1 Gyr, 3Gyr, 
and 10 Gyr. Sources from Table 3 are shown with asterisks. 


appear point-like). Based on the SDSS images and pho¬ 
tometry (see Figure 4, right panel), we also conclude that 
#2 and #32 are likely (bright) foreground stars. The 
SDSS images suggest that #12 and #9 are AGNs (in 
agreement with Figure 4). 

We also attempted to verify the classification using the 
SIMBAD catalog] However, we found that the SIMBAD 
classifications are incorrect for several faint sources. For 
example, we find AGNs (sources #9 and #12, resolved 
with HST ) to be classified as stars. Also, a few stars 
are referred to as QSOs/AGNs in SIMBAD. Therefore, 
we do not rely on the SIMBAD for classifications of the 
X-ray sources in Draco. 

As an aid to help with the classification, we have 
developed an automated machine le arning classification 
pipeline (MUWCLASS; |Brehm et ah|20T4) that we sub¬ 
sequently applied to the X-ray sources m Draco. The 
pipeline uses a learning decision tree algorithm, to per¬ 
form the X-ray source classification (see A1 for details). 


5 http://simbad.u-strasbg.fr/simbad/ 


The automated procedure produced 20 classifications 
with confidence > 70%, including 14 AGN, 5 stars, and 
1 CV (see Table 3). The classifications with lower con¬ 
fidence are indicated by a question mark. Note that the 
calculated confidences for each type do not include the 
uncertainties associated with the X-ray flux determina¬ 
tion which are substantial for faint X-ray sources and 
also affect HRs. They also do not include a possibility 
of assigning false MW counterparts to the X-ray sources 
(confusion), which, however, is expected to be small (see 
Section 3.2). For a number of interesting X-ray sources, 
which could belong to Draco or be XRBs or compact 
objects in our Galaxy, we have carried out a more de¬ 
tailed analysis (see Section 4.1) by extracting the X-ray 
spectra from the original data and performing spectral 
fits (not used in the automated classifications). The lo¬ 
cations of the sources in the diagrams shown in Figures 
3,4 and 5 are mostly consistent with the automated clas¬ 
sifications or, whenever this is not the case, the reasons 
for discrepancy are discussed in Section 4.1 on source- 
by-source basis. 

After we submitted this paper, we became aware of 
two other papers (Manni et al. 2015 and Saeedi et al. 
2015) discussing X-ray source classification in Draco and 
submitted roughly at the same time as our paper. The 
Saeedi et al. (2015) paper presents a more comprehensive 
analysis than Manni et al. (2015) and, bearing in mind 
differences in class definitions, nearly all our classifica¬ 
tions agree with Saeedi et al. (2015), except for source 
15 for which our automated tool preferred AGN (84%) 
classification, while Saeedi et al. (2015) suggest a binary 
star (not XRB) classification (but uncertain). We also 
note that the machine-learning method used in our paper 
is different from those used by Saeedi et al. (2015) and 
Manni et al. (2015) who relied on previously established 
heuristic prescriptions from multiple publications. How¬ 
ever, the very good agreement between our and Saeedi et 
al. (2015) classifications shows that both methods worked 
well for this field and similar results were obtained with 
different classification approaches, adding credibility to 
both independent analyses. 
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TABLE 3 

Properties and automated results for classification X-ray sources in the Draco field 


# 

RA 

Dec 

Sr a 

Flux b 

HR2 C 

HR4 C 

S/N d 

9 e 

i e 

Wl f 

Class 8 (Prob.) 

1 

260.0917 

57.9740 

0.40 

86(2) 

-0.43(1) 

-0.86(1) 

115.4 

12.80 

11.81 

7.73 

GM Dra h 

2 

260.4930 

57.8228 

0.56 

157(16) 

-0.26(2) 

-0.31(6) 

51.3 

11.21 

10.10 

8.25 

Star/CV ? 

3 

260.1406 

58.1417 

0.58 

74(8) 

-0.36(2) 

-0.11(6) 

48.4 

17.34 

16.24 

12.83 

AGN (84%) 

4 

260.1575 

58.0366 

0.49 

21(2) 

-0.24(3) 

-0.68(6) 

43.4 

11.39 

12.19 

7.80 

Star 1 (99%) 

5 

260.2180 

57.9202 

0.47 

18(2) 

-0.20(4) 

0.08(7) 

31.8 

20.29 

19.96 

14.84 

AGN (99%) 

6 

259.8937 

57.9806 

0.47 

15(3) 

-0.33(4) 

-0.1(1) 

26.1 

19.34 

19.02 

13.90 

AGN (99%) 

7 

259.8409 

57.8290 

0.52 

5(2) 

-0.08(8) 

-0.1(2) 

16.1 

22.76 

21.69 

15.98 

? 

8 

259.7567 

58.0081 

0.67 

6(2) 

0.0(1) 

0.2(2) 

14.4 

18.79 

18.68 

15.18 

AGN (97%) 

9 

260.1548 

57.8155 

0.79 

12(3) 

-0.73(6) 

0.5(1) 

14.3 

19.32 

18.28 

13.97 

AGNJ (99%) 

10 

259.9750 

57.9977 

0.71 

l(2) m 

0.91(1) 

0.78(4) 

14.2 

22.97 

21.88 


? 

11 

260.1798 

57.9119 

0.69 

5(1) 

0.5(1) 

0.4(1) 

14 

20.93 

20.22 

15.93 

C V /AGN 

12 

259.7936 

57.8301 

0.75 

10(2) 

0.69(9) 

0.44(9) 

13.4 

18.96 

17.46 

13.87 

AGNJ (88%) 

13 

259.7689 

57.7914 

0.81 

5(2) 

-0.4(1) 

0.3(2) 

12.5 

20.49 

19.12 

15.75 

AGN (99%) 

14 

260.2161 

57.6997 

0.86 

15(5) 

0.1(1) 

0.3(2) 

12.5 

19.85 

19.72 

15.41 

AGN (99%) 

15 

260.0755 

57.8516 

1.65 

2.6(8) 

-0.3(1) 

-0.2(2) 

12.4 

19.48 

16.85 

13.88 

AGN (84%) 

16 

259.8326 

57.9952 

0.85 

9(2) 

0.9(1) 

0.61(8) 

12 

18.80 

17.68 

14.97 

AGN ? 

17 

259.9902 

57.8351 

0.86 

0.1(4) m 

-1.0(l°' 3 ) n 

0-9(iol) n 

10.6 

17.75 

16.19 

13.26 

Cl k 

18 

260.3452 

57.8417 

1.02 

5(2) 

-0.0(1) 

0.1(2) 

10.4 

18.05 

17.81 

14.89 

AGN (84%) 

19 

259.7686 

58.0586 

0.86 

6(3) 

-0.2(1) 

-0.2(3) 

10 




? 

20 

260.4225 

57.8769 

1.11 

9(3) 

-0.3(1) 

0.3(2) 

10 



15.75 

AGN? 

21 

259.9869 

57.7715 

1.17 

5(2) 

0.2(1) 

0.2(2) 

9 




? 

22 

259.7530 

57.8628 

0.99 

2(1) 

0.0(1) 

-0.2(3) 

8.8 

20.27 

19.18 

16.19 

AGN (99%) 

23 

259.9394 

57.8825 

1.05 

0.9(5) 

-0.1(1) 

-0.2(3) 

8.6 

23.20 

21.33 

16.01 

? 

24 

259.8536 

57.7479 

1.14 

2(1) 

0.1(2) 

-0.3(4) 

8.3 




? 

25 

260.1971 

57.8889 

1.17 

1.3(9) 

-0.0(2) 

0.0(4) 

7.9 




? 

26 

259.9040 

57.7832 

1.16 

1.3(9) 

-0.4(2) 

0.1(3) 

7.9 

22.84 

21.35 


AGN? 

27 

260.2515 

57.8649 

1.05 

1(1) 

-0.0(4) 

-0.1(2) 

7.6 

21.45 

21.53 


AGN? 

28 

260.3006 

57.8342 

1.06 

3(2) 

-0.2(2) 

0.4(2) 

7.6 

22.99 

23.18 


? 

29 

260.1186 

57.9400 

1.20 

0.6(5) 

0.4(2) 

0.0(3) 

7.3 




? 

30 

260.4516 

57.9685 

2.00 

2(3) m 

-0.4(2) 

0.1(5) 

7.2 

20.22 

19.81 

16.19 

AGN (99%) 1 

31 

260.2698 

57.8927 

1.21 

2(2) 

-0.2(2) 

-0.1(4) 

7.2 

15.97 

15.03 

13.17 

Star (99%) 

32 

260.1072 

57.8848 

1.16 

0.3(3) 

-0.9(1) 

-i-o(±g- 5 ) n 

7.1 

14.30 

11.30 

10.16 

Star (95%) 

33 

259.9014 

57.8166 

1.18 

0.9(6) 

-0.0(2) 

-0.2(3) 

6.8 




? 

34 

259.6866 

57.8509 

6.27 

6(6) 

-0.5(2) 

0.3(4) 

6.2 




? 

35 

260.1380 

58.1042 

1.69 

1(1) 

0.0(2) 

-0.4(6) 

6.1 




? 


Note. — a la positional uncertainty in arcseconds (from 3XMM-DR5 catalog). ^Observed X-ray fluxes in the 0.5 — 12 keV range 
in units of 10 —14 erg s —1 cm -2 . The value in brackets is the measurement uncertainty in the last digit. c Hardness ratio calculated as 
HR2 = (F( X 2-2 — F 0 .2-i.2)/(Fi. 2-2 + Fo. 2 - 1 . 2 ) and HR4 = (F 2 _7 — Fo.5-2)/(F 2 _7 + Eo. 5 - 2 ), where F X -y is the observed flux in the 
X — Y keV energy band. The value in brackets is the measurement uncertainty in the last digit. d Signal-to-noise ratio. e g and i magnitudes 
front the USDSS DR 9 catalog. ^W1 magnitudes front the WISE catalog. ^Classification confidence (see text). ^GM Dra, var iable star in 


Tycho 2 catalog ( Hpg et al.|2000| ). l A cool K9 type star with X-ray luminosity of 3.7 x 10 29 ergs s 1 (a t a distance of 125 pc; 


al.|2006). 

morpHoIo 


Ammons et 


^Visually resolved in the SDSS images. k Symbiotic star in Draco, see Belczynski et al. (2000). ^Source #30 also shows complex 
morphology (likely, superposition of star and galaxy as can be seen in Figure 2. ""Three sources have the uncertainties of fluxes in the 
individual bands (1—2 keV, 2—4.5 keV, etc.) larger than the fluxes themselves (according to 3XMM-DR5), which results in the large total 
flux uncertainty. n In these cases the uncertainty is strongly asymmetric, because the value of HR hits the boundary, and therefore two 
uncertainties are provided. 


3.2. Cross-correlation and chance superposition 

We cross-correlated all X-ray sources listed in Table 3 
with the USNO-B1, SDSS, 2MASS, and WISE catalogs. 
We find a single optical/NIR/IR counterpart within the 
2" search radiu^] for all but one sourc^] X-ray source. 
The only sources with multiple counterparts within the 
search radius are source #5, which has two USNO-B1 
counterparts, and source #32, which also has two coun¬ 
terparts in SDSS. For those two sources we only consider 
the brightest (which also happens to be the closest) coun¬ 
terpart to be the “true” counterpart below. We have 
verified that all optical/NIR/IR counterparts are within 

6 This happens to be an upper bound on la position uncertain¬ 
ties for 34 out of 35 X-ray sources in the Draco field based on the 
3XMM-DR5 positional parameters (see Table 1). 

7 This is source #34 which is both very faint and strongly off- 


^0"5 of one another. We see no systematic offset be¬ 
tween the X-ray and 2MASS sources and the RMS offset 
between the X-ray and 2MASS/SDSS positions for the 
35 sources listed in Table 3 is ~ 1", i.e. well within the 
distribution of uncertainties for the X-ray source posi¬ 
tions (see Table 3). We note that while the HST images 
reveal multiple counterparts for some of the sources, the 
surveys are not as deep, and, therefore, pick up only 
the brightest objects in the 2" search radius around each 
source. Unfortunately only 5 X-ray sources are located 
within the HST fields and even for them multiband HST 
photometry is lacking. We cannot exclude that some of 
the brighter X-ray sources are associated with the much 
fainter optical sources that can only be seen in the HST 
images. However, these would have to be sources of a 
very rare type rare types (e.g., solitary neutron stars or 


axis. 
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quiescent BH binarie^] in Draco or extended Galactic 
halo) or yet unknown classes. As for the fainter X-ray 
sources, we can hardly say much without deep XMM- 
Newton and HST observations. 

For each survey we calculate the chance superposi¬ 
tion using three different techniques. First, based on 
the average optical/NIR/IR source densities in the field 
(p = 0.0048, 0.01, 0.0004, and 0.001 stars/arcsec 2 in 
the USNO-B1, SDSS, 2MASS, and WISE surveys, re¬ 
spectively) we calculate the probability of finding find¬ 
ing zero field sources in the circle of radius r = 2", 
P = exp(— pirr 2 ). This leads to chance coincidence prob¬ 
abilities of 1 — P — 6%, 12%, 0.5%, and 1% (for USNO- 
Bl, SDSS, 2MASS, and WISE surveys, respectively). For 
the second method we run 1,000 simulations where we 
randomly populate with 35 synthetic X-ray sources the 


optical/NIR /IR field of Draco (similar to e.g., Rangelov 


et al.||2 012). We then count the average number of optU 
cal/N1R/1R counterparts to the synthetic X-ray sources 
within the r = 2" radius. This method results in 0.9, 
1.9, 0.08, and 0.27 (for USNO-B1, SDSS, 2MASS, and 
WISE surveys, respectively) spurious detections of the 
35 sources. The third method we employed was to offset 
each X-ray sources in random direction (within 100" of 
its original position). This led to 1, 1.8, 0.1, and 0.4 (for 
USNO-B1, SDSS, 2MASS, and WISE surveys, respec¬ 
tively) spurious detections in 1,000 simulations. There¬ 
fore, all three methods give very similar results suggest¬ 
ing only 1-2 spurious detections, compared to the 26 X- 
ray sources for which we find optical/NIR/IR counter¬ 
parts (see Table 3). 


4. RESULTS AND DISCUSSION 

Our goals are to identify and/or set constraints on the 
numbers of XRBs that may belong to Draco, set a limit 
on the emission from a possible central intermediate mass 
BH, and look for a possible DM decay line in the diffuse 
X-ray emission. 


4.1. X-ray sources 

Below, we discuss individually several sources that do 
not appear to be AGNs or foreground Galactic stars 
based on the color-color and color-magnitude diagrams, 
HST/SDSS images, and automated classification. For 
the brightest of these sources we performed spectral fits 
with blackbody radiation (bbodyrad), disk blackbody 
(diskbb) and power law (PL) models modified by in¬ 
terstellar absorption (phabs model in XSPEC). In gen¬ 
eral these models provide an adequate description of the 
thermal and non-thermal components of X-ray spectra 
typically observed in the 0.2—lOkeV range for objects 
such as CVs, isolated neutron stars, and LMXBs. The 
best-fit parameters for these sources are given in Table 3. 

Source #10: The optical-IR properties of this source 
exhibit characteristics of an old (few Gyr), evolved star. 
The X-ray spectrum is well fitted with the bbodyrad 
model with a temperature of of 1.11 db 0.08 keV. Nor¬ 
malization of the bbodyrad model corresponds to an 
effective emission region of 0.44 km at 80 kpc. The 
absorption n# = 3.5 x 10 21 cm -2 exceeds the Galac¬ 
tic tih = 2.69 x 10 20 cm -2 in this direction ( Chandra 


Colden toolkit]^]). The source diskbb model gives a sim¬ 
ilar tih and temperature of 2.1 ±0.3keV. If the source is 
in Draco, then the normalization of the diskbb model, 
{{Rin/km)/(D /10 &pc)) 2 cos0, 6 is the angle of the disk 
(assumed face-on), corresponds to r^ n = 0.11 km, which 
seems to be too small for an accretion disk in an XRB. 
The X-ray luminosity, Lx — 8 x 10 34 erg s -1 , is too high 
for a redback or black widow type binary unless it goes 
through an ou tburst episode and sw itches to an accretion 
state (see e.g., Stappers et al. 2014). However, the X-ray 
lightcurve does not show any significant variability dur¬ 
ing all five XMM-Newton observations. Alternatively, it 
is possible that #10 is a non-accreting MSP in a Galac¬ 
tic binary with a late type companion at a distance of 
~ 1 kpc with the X-rays being emitted from the hot po¬ 
lar cap of the NS. However, the large (compared to the 
Galactic) absorbing column would be puzzling in this 
scenario. 

Source #11: The X-ray spectrum requires two com¬ 
ponents, bbody (BB)+PL, as PL or BB models alone do 
not provide acceptable fit to the data. Both the tem¬ 
perature (kT = 0.7 ± 0.1 keV) and the photon index 
(r = 1.2T0.5) are not uncommon for several object types 
(e.g., AGN and CV) but even for the two-component 
model the data show systematic excess below 0.4 keV 
over the best fit, suggesting that even softer component 
is needed. Although the luminosity of 4.5 x 10 32 erg s -1 
at d = 80 kpc is typical of a redback or black widow 
system, the spectrum is unusually soft for these kind of 
binaries. The absorbing column is uncertain and low, 
consistent with the Galactic value. The HST observa¬ 
tions (Figure 2) show a point-like source and an extended 
structure, which may or may not be related to the point 
source. Given the large density of background galaxies 
in the deep HST images, it is possible that we are see¬ 
ing the superposition of a star (either in Draco or in the 
Galaxy) and a background galaxy/AGN. This could ex¬ 
plain why the automated classification is confused (split 
between AGN and CV). We found no variability between 
five XMM-Newton exposures. 

Source #16: This source is strongly (intrinsically) 
absorbed with an odd jump > 8 keV in the X-ray spec¬ 
trum. The source’s optical/NIR properties are consistent 
with those of an evolved low-mass single star, but the 
apparently high absorption in X-rays rules out an active 
corona in a nearby star as a source of X-rays. It is possi¬ 
ble that this source is a symbiotic star or CV in Draco. In 
fact, the X-ray spectrum shows signs of emission line(s) 
from the iron complex (Figure 6). The spectrum also 
shows hints of other lines at lower (1—2keV) en ergies. 


This could b e indicative of a magnetic CV (Fujimoto & 
Ishida|19~97 ). We would like to note that symbiotic stars 


are not a separate class in our current training dataset 
for MUWCLASS and therefore such a system would be 
classified as another object type. This source is classified 
as AGN by the automated algorithm but with low confi¬ 
dence (71%). The X-ray lightcurve does not exhibit any 
significant variability. 

Source #17 : This source has been report e d as a can- 
didate symb iotic star (Belczynski et al. 2000). [Aaron son 


et al. (1982) identified this source as carbon star with 


an unusual SED showing strong emission lines. The au- 


8 such a Swift J1357.2—0933 


9 http://cxc.harvard.edu/toolkit/colden.jsp 
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TABLE 4 

Best-fit spectral results for selected point sources in the Draco field. 


Source # 

Model 

Normalization 8, 

n H h 

kT, Tj n e 

r 

xl 

«/ f 




(10 22 cm -2 ) 

(keV) 




#10 

PL 

(1.534=0.48) x 10 -5 

1.35T0.27 


1.72T0.23 

0.85 

171 


diskbb 

(1.92±1.15)xl0 -4 

0.83T0.16 

2.09T0.35 


0.85 

171 


bbodyrad 

(3-354=0.89) x 10 -3 

0.35T0.14 

1.11T0.08 


0.88 

171 

#n 

bbodyrad+PL 

(5.644=2.94) xlO -3 c 
(2.20±1.35)xl0 -6 d 

< 0.027 

0.75T0.11 

1.24T0.48 

1.05 

171 

#16 

PL 

(1.854=0.50) xlO -5 

0.92T0.17 


1.82T0.21 

1.16 

60 

#19 

PL 

(1-014=0.16) xlO -5 

0.06T0.04 


2.91T0.38 

0.77 

43 


bbody 

(3-324=0.28) xlO -7 

< 0.027 

0.18T0.01 


1.06 

43 


bbodyrad 

2.21±0.63 

< 0.027 

0.18T0.001 


1.06 

43 

#23 

PL 

(3.47±1.98) x 10 -6 

0.23T0.18 


3.6T1.3 

0.89 

75 


bremss 

(5.634=5.56) x 10 -6 

< 0.027 

0.83T0.55 


0.88 

75 

#25 

diskbb 

(2.124=0.14) x 10 -3 

< 0.027 

0.64T0.11 


0.83 

197 


a Normalization of the bbodyrad model is R^ m /.D 2 0 , where R^ m is the source radius in km, Dio is the distance to the source in units of 
lOkpc and the normalization of the diskbb model given as ((Ri n /km)/(D/10 kpc)) 2 cos$, where R* n is the is ”an apparent” inner disk 
radius. 

b Galactic extinction is assumed as 2.7xl0 20 cm -2 for each source. 
c Normalization of the bbody component. 
d Normalization of the PL component. 

e kT; temperature keV from bbody and bbodyrad models. Tj n ; temperature at inner disk radius (keV) from diskbb model. 
f number of degrees of freedom. 


thors suggest that the optical colors can be explained by 
a symbiotic binary (red giant with hot main sequence 
companion). Aaronson et al. (1982) also claim that #17 
is in Draco based on radial velocity measurements. How¬ 
ever, the observed X-ray flux implies a high X-ray lumi¬ 
nosity of Lx = 1.2 x 10 33 erg s -1 at the distance of Draco, 
which can hardly be produced in a non-degenerate binary 
with red giant and late-type main-sequence star. There¬ 
fore, we consider this source to be a good candidate for 
a quiescent XRB in Draco. The automated classification 
did not produce a confident result. The confusion could 
be due to the unusual SED noticed by [Aaronson et al. 


(1982) who also suggested that the star can be in a bi- 


nary with a (pre-)degenerate object or it can be a star 
caught very early in the process of ejecting the outer lay¬ 
ers and forming a planetary nebula. If the object does 
not fit any of the predefined object classes the automated 
algorithm is expected to be confused. 

Source #19: Given its X-ray flux and lack of op¬ 
tical counterpart (arguing against binary nature), it is 
unlikely that this source would belong to Draco (X-ray 


dence of even further hardening beyond 5—6 keV in pn. 
PL and bremss provide a reasonable fit to the data. Our 
automated algorithm does not provide a confident classi¬ 
fication, likely due to incomplete MW parameters (a very 
faint counterpart found in SDSS, but not in 2MASS). 

Note that nine sources lack counterparts in SDSS. 
These could be AGNs that are too faint to be detected 
in SDSS , or unusual objects that produce little emis¬ 
sion in the optical, such as isolated NSs or quiescent BH 
binaries in Draco. Of those nine, only source #25 was 
imaged by HST (Figure 2). Two point-like objects are 
seen within the r = 2" circle, one of which could be the 
counterpart of X-ray source. If this is the case, it could be 
LMXB/CV in Draco (based on the optical brightness). 
Unfortunately, no multi-band photometry exist for this 
field, which prevents us from classifying the possible op¬ 
tical counterpart of the source. 

After careful examination of all X-ray sources we find 4 
potentially interesting sources (#10, #16, #17 and #25) 
that could belong to Draco and be potential LMXBs, 
CYs and/or symbiotic candidates. We note in passing 


is too high for isolated NSs), so it could be an old re- 

et al. 2005 

) and Fornax to have 2—3 “field” LMXBs (i.e., 

cycled or non-recycled pulsar in a non-accreting binary 

not part 01 

■ Fornax globular clusters, which have another 

in our Galaxy. However, the source is imaged strongly 

2 LMXBs; 

Nucita et al. 2013 

). Given the small number 


positional uncertainty, even though source #19 has one 
of the smallest positional uncertainties in Table 3 (pos¬ 
sibly underestimated by the automated source detection 
tool in 3XMM). There is an SDSS source within 3" with 
u = 20.6, g = 20.4, r = 20.3, i = 20.3, and 2 = 20.3. If 
this SDSS source is a counterpart of the X-ray source, 
it would exclude the possibility of #19 being nearby 
isolated pulsar/NS, which otherwise would be possible 
based on X-ray spectrum and lack of optical/IR coun¬ 
terparts. The lack of optical/IR data prevents our algo¬ 
rithm from confidently classifying this source. 

Source #23: The source is faint with a soft X-ray 
spectrum, too soft for typical AGN. There is some evi- 


(except, perhaps, source #16 and #17), our findings are 
roughly on par with those for Sculptor and Fornax (both 
of which are Milky Way satellites and at roughly the same 
distance as Draco). However, we caution the reader that 
the respective analyses use different techniques and rely 
on different MW data and therefore a direct comparison 
is speculative at best. 


Walker et al. (2015) recently presented optical spec¬ 
troscopy for stars m the Draco field. They fit stellar 
models to the all sources and produced a dataset with 
the derived properties (line-off-sight velocity, VJ os , ef¬ 
fective temperature, T e ff, surface gravity, g, metallicity, 
[Fe/H], etc.). We cross-match the spectroscopic dataset 























with the positions of the 35 X-ray sources, and found 5 
matches - sources #9, #16, #17, #22, and #30. The 
optical spectroscopy for sources #9, #22, and #30 has 
very low signal-to-noise (S/N < 3), which resulted in 
very uncertain fit par ameters and, hence, no classifica¬ 
tions are provided in Walker et al.| (2015]). These low 
signal-to-noise sources were classified (by our algorithm) 
as AGN (in agreement with Figure 4). Sources #16 and 
#17 have Vi os = —295 db 1, which is consistent with Vi os 
of Draco, and therefore, suggests the association of the 
two sources with the dSph. Spectroscopy suggests that 
these two sources are evolved, low-mass stars, which is 
consistent with CV/LMXB scenarios discussed above. 


4.2. Limits on central IMBH mass 

It has been speculated that dSphs can host an inter¬ 
mediate mass BH s (IMBHs; |Nucita et al. 2013 Mac- 


carone et al.||2005| and references therein). We estimate 


the X-ray luminosity for a putative IMBH in Draco due 
to Bondi-Hoyle accretion as 


L x ^ 3.4 x 10 37 e 


Mbh 


10 3 M 0 


3 ) (lOkms-O 

( n . ~ _ 3 ) erg s' 
V0.1 cm 6 J 


-3 


(i) 


where e is the radiative efficiency, and n is the hydrogen 
number dens ity in D raco . Based o n observational results 
( Mateo|1998 Grcevich fc Putman||2009 ), we estimate an 
upper limit of n < 0.02 cm - " T Jand sound speed c s « 
10 km s _ 1 for dSphs, whic h should be of the order of the 
stellar veloc ity dispersion (Grebel et al.|2003 Hargreaves 
et al. 1996). Therefore the corresponding flux at the 
distance of Draco is 



Mbh ( - ”V /2 < 3.4 x 10 3 ( by ) 

\0.02 cm -3 / V10- 15 cgs J 

3/2 


1/2 


( e \-V 2 / ^ 

(lFs) (lO km s-0 M ®' 

(3) 

Theoretical models a nd observational constraints (e.g., 


Gebhardt et al. 2000) show that the central BH mass 
would depend on the velocity dispersion as 


o / a \ 3 - 75 

M bh = 1.6x10 ( 10 km ^=t) M ®’ ( 4 ) 

based on the Draco velocity dispersion measurement by 
Hargreaves et al. (1996). This mass estimate is compa¬ 
rable with the limit calculated above. We note, however, 
that the mass estimate is based on the use of the upper 
limit on the hydrogen density and a rather optimistic ef¬ 
ficiency of 10 -5 . More conservative values for the density 
and (especially) the efficiency (where the uncertainty is 
the largest) could easily reduce the mass limit by several 
orders of magnitude. 


4.3. DM decay feature at 3.5 keV? 

We also searched for the ^3.5 keV line that has been 


previously report ed in the Perseus galaxy clu ster (Bo- 


yarsky et al.|2014) and M31 (Bulbul et al. 2014). Figure 7 
shows the X-ray spectrum of the diffuse EPlG-pn emis¬ 
sion from 0.5 to 10.0 keV. No significant spectral features 
consistent with previously reported line at ^ 3.5 keV is 
seen. 


5. SUMMARY 

We have analyzed 35 (bright) X-ray sources detected 
by XMM-Newton in the Draco dwarf spheroidal galaxy. 
Along with traditional techniques of classification such as 
the use of color-color and color-magnitude diagrams and 
hardness ratios, we have utilized an automated machine¬ 
learning approach to identify X-rays sources in Draco. 
These were complimented by the analysis of HST /SDSS 
images and X-ray spectra for several interesting X-ray 
sources. Our main results are: 


According to the 3XMM catalog there are four X-ray 
sources with S/N ~ 3 — 5 within ^ 3' of the center of 
Draco (FI.A. = 17 /l 20 m 13.2 s and Decl. = +57°54 / 55'/3; 
Rave et ah] [2003). The measured flux for the nearest 
source is Fx ~ 4.5 x 10 -15 ergs -1 cm -2 (the other 
th ree sources have similar fluxes and S/N). According 
to Beskin & Karpov (2005) the expected radiative effi¬ 


ciency for synchrotron radiation (presumably the domi¬ 
nate radiation mechanism for a slowly accreting BH) is 
e ^ 10 -5 . Assuming that one of these X-ray sources is 
the putative IMBH, we can then estimate its mass. How¬ 
ever, it is possible that none of these faint X-ray sources 
represent the putative IMBH. Therefore, we also can de¬ 
termine the lower limit on the IMBH flux by measuring 
the diffuse flux fluctuations (^ 10 -15 erg s -1 cm -2 at 
3 a level) in the local background in the combined EPIC 
image. This leads to the following constraint 


• The classification of X-ray sources in Draco re¬ 
sulted in 12 AGN and 3 foreground stars. We 
also identified 4 X-ray sources (potential quiescent 
LMXBs or CVs) that could belong to Draco. For 
two of them (#16 and #17) the associations are 
supported by the line-of-sight veloci ty measure¬ 
ments reported in Walker et al.l (2015). 


• The upper limit on the mass of the IMBH that we 
obtained from the X-ray data analysis is similar 
to the predictions based on the velocity dispersion 
correlation. The current data are consistent with 
the none existence of a central black hole, however 
a deeper observation may result in such a detection. 




We do not find any significant evidence of 3.5 keV 
emission line which could be a signature of decaying 


DM. 
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Fig. 6.— Best fit model spectra and their residuals for sources #10, #11, #16, #19, #23, #25 which are not classified as AGN or 
foreground galactic stars by our automated classification pipeline. MOS1, MOS2, and EPIN-PN data points and their respective best 
model fits are shown in black, red and green. 












































































































































































































































































































































































































































































10 



Fig. 7. — The combined (all five observations) EPIC-pn spectrum 
of the diffuse X-ray emission. The background spectrum for Draco 
is shown in red in the 0.5-10.0 keV energy range. The filter wheel 
closed spectrum for the XMM-Newton EPIC-pn camera is shown in 
black. The broad line features are known to be due to instrumental 
lines; A1 K a at 1.49 keV and Cu florescent at 8 keV. No evidence 
is found for a spectral feature at ~ 3.5 keV. 


in Draco. Currently, we have 9 sources that lack opti- 
cal/NIR detections and may represent solitary NSs or 
quiescent BH binaries in Draco or extended Galactic 
halo. Combined high-resolution multi-band optical ob¬ 
servations, sensitive optical spectroscopy, and a deep 
XMM-Newton observation can provide a much better 
view into the population of X-ray sources in Draco, mass 
of the putative IMBH, and the origin and evolution of 
dSphs. 
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a summer visitors program and The Science Academy 
(Bilim Akademisi, Turkey) under the BAGEP program. 
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ter, which is operated by the Smithsonian Astrophysical 
Observatory for and on behalf of the National Aeronau¬ 
tics Space Administration under contract NAS8-03060. 


Deep multi-band optical observations are required for 
a more detailed study of the X-ray source population 

APPENDIX 

Al. Machine-learning approach to X-ray source classification 


Our automated classification tool, MUWCLASS (Brehm et al. 2014; Hare et al. in prep.), uses C5.CPj which is a 
supervised decision tree learning algorithm requiring a training dataset with already classified sou rces (seebelo w). C5.0 
is an updated version of C4.5 (Quinlan|1993), which is based on the ID3 algorithm created by Quinlan (fl986). We ran 
C5.0 classification algorithm with the default parameters as described in https://www.rnlequ est.com/see5-in fo.html. 


2003) 


We have adopted a Laplace prescription for the estimation of the classification confidences (Chawla et al. 
each leaf, P = (TP + 1)/(TP + FP + C), where TP, PP, and C are true positives, false positives, ana number 
of classes, respectively, for the leaf where the source in question has ended up. Currently, MUWCLASS does not 
take into account the uncertainties of the MW parameters while performing the classification. Our machine-learning 
classification approach in this paper closely follows that of Lo et al. (2014) and Farrell et al. (2015) except that we 
do not use the X-ray variability parameters since none of the Draco sources appears to be variable at the statistically 
significant level. 

The training dataset is used to evaluate the parameters of objects from known classes and build the decision tree. 
The decision tree, calculated from the training data set, is applied to a set of unclassified objects. We have invested 
significant effort to investigate the catalogs and related literature to create a training dataset of ~ 8,000 confidently 
classified X-ray sources detected by either Chandra or XMM-Newton (Hare et al. in prep.). This dataset includes 
dif ferent num b ers of sources from various classes. To compensate for the imbalance we use the SMOTE technique 
(jChawla et al.||2011|) to create an expanded balanced dataset. The dataset we used here has 9 predefined objects 
classes: (1) mam sequence stars (General Catalog of Variable Stars; Samus et al. 2009), (2) young stellar objects 
(Chandra Orion Ultradeep Point Source Catalog and PAN-Carina; Getman et al. 2005, Povich et al. 2011), (3) AGNs 
(Veron Catalog of Quasars & AGN; Veron-Cetty & Veron 2010), (4) LMXBs (Low-Mass X-Ray Binary Catalog, 2007; 
Liu et al. 2007), (5) HMXBs (Catalog of High-Mass X-Ray Binaries in the Galaxy; Liu et al. 2006), (6) cataclysmic 
variables (CVs; Cataclysmic Variables Catalog, 2006, Downes et al. 2001), (7) isolated neutron stars (NSs; ATNF 
Pulsar Catalog; Manchester et al. 2005), (8) binary non-accreting NS (ATNF Pulsar Catalog), and (9) Wolf-Rayet 
stars (The Vllth Catalog of Galactic Wolf-Rayet Stars; van der Hucht 2001). We cross-match all objects from the 
training dataset with MW data and extract the following parameters: X-ray fluxes in 4 bands from the 3XMM-DR5 
catalog (the same as those defined in Section 3.2), optical ugriz magnitudes form SDSS, NIR jhk magnitudes from the 
Two Micron All-Sky Survey (2MASS), and IR Wl, W2 and W3 magnitudes from Wide-field Infrared Survey Explorer 
(WISE). 

Figure 8 shows the cross-validation matrix for the SMOTEed training dataset which has the true object class on 
X-axis and the inferred class on Y-axis. The more diagonal is the matrix, the better is the performance of the 
classification algorithm. We have also verified the classification performance by dividing the dataset in two parts and 
using the first part for the classification tree construction and second half for the validation. The dataset is divided 


10 http://www.rulequest.com 
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Fig. 8. — The cross-validation matrix of the training dataset (see Al). Original classes are given on the X-axis, and the classification 
outcome is on the Y-axis. NS_BIN refers to binary non-accreting NS. 


into parts before the SMOTE procedure is applied to the part that is used for training. We typically find that ~ 93% 
of the sources are classified correctly, with the AGN and stars being found more confidently than the other sources. 
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